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Abstract.  We  consider  the  n(n— 1)/2  uniquely  defined  positive 
intervals  among  the  first  n  <  106  prime  numbers  as  a  probe  of 
the  global  nature  of  the  sequence  of  primes.  A  statistically  strong 
periodicity  is  identified  in  the  counting  function  giving  the  total 
number  of  intervals  of  a  certain  size.  The  nature  of  the  periodic 
signature  implies  that  the  sequences  of  intervals  spanning  fixed 
numbers  of  gaps  repeat  quasi-cyclically.  From  the  distribution  of 
intervals  we  extract  also  the  characteristic  period  of  the  repetition, 
which  increases  with  n  in  a  step-wise  manner  between  consecutive 
primorial  numbers  and  coincides  with  the  most  commonly  occuring 
interval.  The  relationship  between  the  most  common  interval  and 
the  primorial  numbers  is  noteworthy  independently  of  the  periodic 
behaviors. 


1.  Introduction 


The  sequence  p  =  2,  3, . . .  of  prime  numbers  exhibits  a  perplexing  hy¬ 
brid  of  locally  unpredictable  but  globally  regular  behaviors  [10].  The 
Prime  Number  Theorem  (PNT)  represents  the  most  significant  global 
property  of  the  sequence.  The  PNT  may  be  stated  concisely  in  terms  of 
the  prime  counting  function  n (y),  defined  to  give  the  number  of  primes 
no  greater  than  y,  as 
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It  follows  from  Eq.(l)  and  n  =  tt (pn)  that 

(2)  pn  ~  nlog(nlog(n)) , 

where  pn  is  the  n-th  largest  prime  and  denotes  equality  in  the  limit 
as  n  diverges.  Throughout  the  remainder  of  this  work  n  is  implicitly 
very  large. 

The  global  behavior  identified  in  the  PNT  may  be  interpreted  as  an 
emergent,  aggregate  property  of  the  sequence  gi,g2, . . .  of  prime  gaps, 
where 

(3)  gn  Pn+\  Pn 

is  the  n-th  prime  gap.  From  Eqs.(l)  and  (3)  we  have 

(4)  gn  ~  —  ~  log(pn), 

n 

which  implies  that  the  average  gap  in  the  vicinity  of  a  given  pn  is 
approximately  log(pn).  The  local  distribution  of  gaps  with  respect 
to  the  average  is  irregular  and  similar  in  nature  to  a  probabilistically 
determined  distribution  [3],  [6],  [5].  An  effective  measure  of  the  global 
distribution  of  gaps  is  the  counting  function 

(5)  Gp(n,  g)  =  #  {pm  <  pn  :  gm  =  g} 

giving  the  number  of  gaps  of  a  given  size  g  among  the  Erst  n  primes. 
G p  =  Gp(n,g)  naturally  exhibits  a  period-6  oscillation  with  respect  to 
g  for  a  given  n.  A  representative  illustration  of  the  periodicity  is  shown 
in  Fig.(l),  which  contains  a  plot  of  Gp(106,g)  for  all  even  g  from  2  to 
44.  Similar  periodicity  is  observed  in  association  with  the  differences 
between  consecutive  gaps  [4], 

It  is  historically  noteworthy  and  epistemologically  significant  that  the 
PNT  evolved  from  empirically  motivated  conjectures.  In  particular,  no 
later  than  1798  Legendre  surmised  from  tabulated  data  a  logarithmic 
property  that  is  asymptotically  equivalent  to  Eq.(l).  Modern  computa¬ 
tional  capabilities  have  facilitated  dramatically  such  empirical  studies 
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Figure  1. 

of  the  properties  of  the  primes.  See.  for  instance,  references  [1],  [2], 
[4], [7],  [8]  and  [9].  Although  it  is  not  generally  recognized  as  a  legit¬ 
imate  field  within  the  mathematical  community,  the  empirical  study 
of  the  nature  of  numbers  has  been  perhaps  as  important  for  the  ad¬ 
vancement  of  number  theory  as  experimental  physics  has  been  for  the 
advancement  of  physical  theory. 

The  purpose  of  the  work  is  to  demonstrate  empirically  a  new  global 
property  of  the  primes  whose  nature  is  both  self-similar  and  periodic. 
Rather  than  consider  only  the  prime  gaps,  we  consider  here  the  entire 
inventory  of  positive  intervals  among  the  first  n  primes.  In  the  fol¬ 
lowing  Section  we  present  a  simple  generalized  formalism  for  analyzing 
the  distribution  of  intervals  in  an  arbitrary  increasing  sequence  of  real 
numbers.  In  Section  3  the  formalism  is  applied  to  the  intervals  among 
the  first  n  primes,  for  n  as  large  as  106.  The  distribution  exhibits  a 
strong,  quantifiable  periodicity  which  is  identified  as  the  signature  of 
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cyclical  repetition  in  the  sequence  of  intervals  spanning  a  given  num¬ 
ber  of  prime  gaps.  Moreover  we  find  that  the  characteristic  period  of 
the  repetition  is  nearly  always  the  most  common  interval  size,  both  of 
which  increase  with  n  in  a  step-wise  manner  and  are  always  primorial 
numbers  (he.  2,  6,  30, 210, 2310, . . . ).  Section  4  contains  a  summary  of 
the  primary  conclusions  followed  by  several  related  conjectures. 

2.  Positive  intervals  within  an  increasing  real  sequence 

Let  x  =  Xi,  X2,  ■  ■  ■  be  an  increasing  sequence  of  real  numbers  which  may 
be  either  finite  or  infinitely  long.  Throughout  the  following  every  bold¬ 
faced  lower-case  Roman  letter  represents  a  particular  type  of  increasing 
real  sequence  whose  most  general  form  is  x.  Among  the  first  n  numbers 
in  x  there  are  n2  uniquely  defined  intervals  which  may  be  organized 
naturally  in  the  form  of  a  matrix.  For  any  given  x  and  n  we  therefore 
define  Ax>n  to  be  an  n  x  n  matrix  with  elements  Xk  —  Xj,  where  k 
and  j  designate  respectively  the  column  and  row.  The  elements  along 
the  primary  diagonal  j  =  k  vanish  trivially,  leaving  only  n(n  —  1) 
uniquely  defined  non-vanishing  intervals.  Furthermore  because  Ax  ,n  is 
antisymmetric  a  complete  description  of  the  intervals  requires  only  the 
n(n  —  l)/2  positive  elements  populating  the  upper- diagonal  region  of 
AXjn.  Let  the  positive  intervals  be  defined  separately  as 

(6)  d-j.k  =  xk  -  xj  ;  Vj  <  k. 

In  the  interest  of  formality  let  DXi„  be  the  set  of  all  n(n  — 1)/2  uniquely 
defined  dj^  among  the  first  n  numbers  in  a  given  x. 

Throughout  the  following  l  is  reserved  for  the  number  k  —  j  of  gaps 
spanned  by  some  dj ^  £  T>^n  and  is  therefore  never  greater  than  n  —  1. 
For  a  fixed  l  all  intervals  djj+i  lie  on  the  same  diagonal  in  Axn,  to 
which  the  term  7-th  diagonal’  refers  henceforth.  The  term  ‘gap’  is 
used  here  in  reference  to  an  interval  between  consecutive  numbers  in 
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any  given  x,  hence 

(7)  Qj  =  djj+i. 

It  is  important  to  note  that  the  elements  of  a  given  Dxn  are  uniquely 
dehned  but  are  not  necessarily  uniquely  valued.  In  order  to  characterize 
the  distribution  of  positive  intervals  of  various  sizes  among  the  first  n 
numbers  in  a  given  x  we  define  the  counting  function 

(8)  Dx(n,  6)  =  #  {dj,k  G  £>x,n  :  djjk  =  5}  • 

It  is  also  useful  to  define 

(9)  D^(n,  6)  =  #  {djj+i  G  Vxn  :  djJ+i  =  5} 

to  give  the  total  number  of  intervals  of  size  6  which  span  exactly  l  gaps, 
hence 

n—  1 

(10)  Dx(n,6)  =  J2Di](.n,6). 

1=1 

In  order  to  represent  the  most  numerous  interval  size  among  the  first 
n  numbers  n  x  let  Ax  =  Ax(n)  be  dehned  such  that 

(11)  Djn,  Ax)  >  Dx(n,S) 

for  all  5.  In  any  case  where  Dx(n,8 )  has  the  same  maximal  value  for 
more  than  one  different  5  then  let  A x(n)  be  dehned  as  A x(n  —  1). 

The  following  generalized  examples  are  presented  in  order  to  establish 
the  basic  properties  of  Dx  =  Dx(;n,5)  relevant  to  the  present  investiga¬ 
tion.  In  so  doing  it  is  useful  to  specify  an  ordered  set  g  =  {71, 72, , . .  7^} 
containing  the  gaps  that  may  be  found  in  a  particular  sequence.  The 
variable  v  is  reserved  for  the  number  of  elements  in  g.  It  is  also  conve¬ 
nient  to  define 

V 

r  =  ^2 1 m 

m—  1 


(12) 
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r 

v' 

such  that  7  is  the  mean  value  among  all  elements  in  a  given  g.  If  v>  1 
then  the  manner  in  which  the  gaps  are  selected  from  g  determines 
the  nature  of  the  sequence  and  is  specified  accordingly.  Note  that  the 
elements  of  g  are  not  necessarily  listed  in  ascending  order,  nor  are  they 
necessarily  all  uniquely  valued.  If  n  3>  v  and  if  the  elements  of  g  occur 
in  the  sequence  with  approximately  uniform  frequency  then 

(14)  dhn  (n  -  l)y. 


2.1.  Uniformly  spaced  sequence.  Consider  first  the  rudimentary 
case  of  a  sequence  u  =  ui,U2,  ■  ■  ■  of  uniformly  spaced  numbers.  By 
definition  the  associated  set  of  available  gaps  is  simply  g  =  {7}  for 
some  7  =  71.  We  therefore  have 


(15) 


dj.t  =  (k-j)  7. 


All  elements  spanning  l  gaps  are  consequently  identical.  The  total 
number  of  intervals  of  a  given  size  in  T>un  is  simply  the  length  n  —  l  of 
the  diagonal  in  AUjn  on  which  they  are  found,  hence 


(16)  Du{n,5) 


n  —  5/7  for  5  =  7,  27, . . . ,  (n  —  1)7 
0  otherwise 


Note  that  the  non-vanishing  values  of  Du  vary  linearly  with  respect  to 

5. 


2.2.  Sequence  with  randomly  ordered  gaps.  Next  let  us  consider 
a  sequence  r  =  77, 77, . . .  in  which  each  gap  is  selected  in  an  effectively 
random  manner  from  some  g  containing  v  >  1  elements,  each  of  which 
is  equally  probable.  An  interval  dj ^  spanning  more  than  v  consecutive 
gaps  in  r  is,  on  average,  given  by 


(17) 


dj.k  &(k-  j) 7 
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analogously  to  Eq.(15).  If  n  3>  v  then  the  indices  j  associated  with 
intervals  dhk  equal  to  a  particular  5  >  T  are  found  with  uniform  prob¬ 
ability  over  the  range  [1  ,  n  — h/7],  where  5/y  is  the  average  number  of 
gaps  spanned  by  an  interval  equal  to  5.  Moreover,  each  allowed  value 
of  djyk  is  equally  likely  to  be  found  in  association  with  a  given  j.  It 
follows  that  the  average  density  Dr/(n—5/ 7)  is  approximately  constant 
among  all  5  for  a  given  r  and  n>E  We  therefore  have 

(18)  Dr(n,5)  oc  n  —  5/^f, 

where  oc  signifies  global  proportionality.  Equation  (18)  implies  that 
the  non-vanishing  values  of  Dr  for  6  >  T  are  correlated  to  a  single  line, 
approximating  the  precise  linearity  in  Eq.(16). 

As  an  example  let  r'  be  a  particular  sequence  of  the  form  r  with  g  = 
{2,4,6}.  Fig. (2)  is  a  plot  of  all  non-vanishing  Hr/(103, 6).  The  largest 
interval  is  3998,  which  is  expected  from  Eq.(14).  The  plotted  points 
are  strongly  clustered  about  the  line  500  —  <5/8,  with  the  associated 
correlation  coefficient  being  approximately  0.997.  Although  the  details 
of  the  linear  trend  are  immaterial  to  the  present  investigation  it  is 
worthwhile  to  mention  that  the  constant  of  proportionality  in  Eq.(18) 
is  approximately  01/7  in  any  case  where  g  is  of  the  form 

(19)  =  {71,  (7i  +w),  (71  + 2a;),  . . . ,  (7i  +  (i/-l)a;)}. 

In  the  present  example  we  have  u  =  2  and  7  =  4,  hence  Dri  is  broadly 
correlated  to  the  line  n/2  —  5/8. 

2.3.  Sequence  with  cyclically  repeating  gaps.  Finally  let  us  con¬ 
sider  a  third  generalized  sequence  c  =  ci,C2,...  in  which  the  gaps 
repeat  cyclically  such  that 

(20)  9j  =  7n  ;  k  =  1  +  mod(j  — 1,  v) 

for  a  given  g  with  v>  1  and  for  all  j  =  1, . . .  ,  n  —  1.  We  therefore  have 

(21)  9j  = 
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FIGURE  2.  Distribution  of  intervals  in  a  sequence  r'  with 
gaps  randomly  selected  from  g  =  {2,4,6}. 

which  implies  the  more  general  periodicity 

(22)  dj+v,kiv 

Because  every  interval  spanning  is  gaps  is  equal  to  T  we  also  have 

(23)  dj^+v  —  djtk  +  r . 

Every  element  on  the  bis- th  diagonal  in  Ac>n  is  equal  to  bT,  for  any 
positive  integer  b  no  greater  than  (n  —  1  )/is,  and  thus 

(24)  Dc(n,  bT)  —  n  —  bis. 

Furthermore  for  mod(/,  is)  ^  0  and  n  —  l>  is  Eq.(22)  implies  that  the 
sequence  d2,2+i,  ■  ■  ■ ,  dn-ijTl  of  all  intervals  spanning  l  gaps,  i.e.  the 

/-th  diagonal  in  Ac  n,  repeats  after  every  is  terms. 

Because  every  interval  spanning  is  gaps  in  c  is  equal  to  T,  the  most 
common  interval  size  is  also  T  for  sufficiently  large  n.  We  may  prove 
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that  property  by  considering  an  upper  bound  on  Dc  for  5  not  equal 
to  an  integer  multiple  of  T.  For  convenience  let  ac\n,S)  be  the  total 
number  of  intervals  equal  to  a  given  5  in  a  single  cycle  of  v  consecutive 
rows  on  the  l-th  diagonal.  If  mod(<5,  r)^0  then  ai l\n,5)  could  be  no 
greater  than  v  —  1.  The  total  number  of  uniquely  defined  dj ^  ^5  on 
the  l-th  diagonal  could  be  therefore  no  fewer  than  the  number  of  full 
cycles  on  the  diagonal,  hence 

(25)  D<‘>(n,S)  <n-l-  (n~'~ 

Since  no  two  intervals  of  the  same  size  could  occur  on  the  same  row  in 
any  Axn,  the  sum  of  <7c\n,  6)  over  all  /  could  be  no  greater  than  v  —  1. 
Consequently  the  upper  bound  on  Dc  occurs  when  <jc\n,5)  —  u—  1 
and  when  all  intervals  equal  to  5  span  just  one  gap.  From  Eq.(25)  we 
therefore  have 

(26)  Dc(n,S )  <  (n  —  1)  ^ ^  ;  V  mod(5,  T)  ^  0. 

It  follows  from  Eqs.  (24)  and  (26)  that 

(27)  Ac(n)  =  T;  V  n  >  u2  —  u+l. 

Note  that  Eq.(27)  represents  the  most  general  and  thus  the  strongest 
possible  requirement.  In  a  case  where  v  is  not  significantly  greater  than 
the  number  of  unique  values  among  the  elements  of  g.  the  requirement 
for  Ac(n)  =  r  is  typically  just  n  >  v. 


The  period-T  repetition  of  intervals  in  c  also  effects  periodicity  in  the 
A-dependence  of  Dc.  With  k—j+l  Eq.(23)  implies 

(28)  <T%>{n,S)=<T<l*\n,S  +  r) 


and  thus  also 
(29) 


D^inJ+T)  _  DP(n,l) 
n  —  l  —  v  n  —  l 


provided  that  there  are  at  least  u  intervals  spanning  l  +  u  gaps  among 
the  first  n  numbers  in  c,  i.e.  for  n  —  l  —  v  >v.  With  both  sides 
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multiplied  by  n  —  l  —  v ,  and  in  summation  over  all  allowed  /,  Eq.(29) 
becomes 


(30) 


n—2is 

D(c+",(n,6+ r) 

1=1 


n—2is 

Dc(n,S)  -  v  Y 
1=1 


n-l  ) 


If  n—2v  is  at  least  as  large  as  <5/7  then  the  difference  between  the  left 
side  of  Eq.(30)  and  Dc(n,5  +  T)  is  negligible  on  average.  For  all  §  no 
greater  than  (n  —  2z/)y  ~  cn  —  2T  we  therefore  have 


n-2v  /  ^(0/  n\ 

(31)  Dc(n,S+ T)  «  Dc(n,  6)  ,-  v  ^  j- 

Note  that  Eq.(31)  is  valid  for  all  but  a  comparatively  small  range  of 
5  whenever  n  3>  z/.  Furthermore  in  the  first-order  approximation  for 
l<Cn,  which  corresponds  to  5<Ccn,  Eq.(31)  becomes 

(32)  Dc(n,5+ r)«  (l-^c(n,i). 

Equation  (31)  implies  that  Dc  varies  with  S  in  a  cyclical  manner  with 
period  T  for  all  5  spanning  on  average  no  more  than  n  —  2v  gaps.  For 
<5<Ccn  the  correlation  reduces  to  the  simple  proportionality  in  Eq.(32), 
such  that  the  basic  form  of  Dc  over  the  range  (0,T]  is  repeated  in 
progressively  smaller  scale  over  subsequent  periods  of  S. 


As  an  example  let  c'  be  a  particular  sequence  of  the  form  c  with  q  = 
{2,4,6},  hence 
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Figure  (3)  is  a  plot  of  non- vanishing  values  of  Hc/(102,  5).  In  accordance 
with  Eq.(27)  the  most  common  interval  among  the  first  100  numbers 
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Figure  3.  Distribution  of  intervals  in  a  sequence  c'  with 
cyclically  repeating  gaps  2, 4, 6,  2, 4,  6, ... . 

in  c'  is  T  =  12,  with  Dc/(102, 12)  =  97.  The  distinctive  fanned  pattern 
in  Fig.  (3)  is  a  signature  of  the  periodicity  in  Eq.(32).  The  periodic 
structure  for  l  <C  n  may  be  seen  clearly  in  Fig.  (4),  which  is  a  zoorned- 
in  view  of  Fig. (3)  for  6  <  36.  Note  how  the  non- vanishing  values  for 
6  G  (0, 12]  are  repeated  in  proportion  over  the  range  (12,24],  etc. 

3.  Positive  intervals  within  the  sequence  of  primes 

Let  ns  now  employ  the  basic  measures  developed  in  the  preceding  Sec¬ 
tion  to  analyze  the  distribution  of  positive  intervals  among  the  first  n 
primes.  Because  of  the  solitary  odd  prime  gap  g\  =  1  there  are  n  —  1 
odd  ditk  €  Dp ,n,  each  of  which  is  uniquely  valued.  The  anomalous  odd 
intervals  are  not  important  in  the  present  investigation  and  are  ignored 


SCOTT  FUNKHOUSER,  CARLA  COTWRIGHT,  DIPENDRA  SENGUPTA,  AND  BRYAN  WILLIAMS 


0  5  10  15  20  25  30  35 

s 


Figure  4. 

throughout  the  following.  Although  it  is  not  mathematically  meaning¬ 
ful  to  construct  p  in  terms  of  some  g,  for  notational  convenience  let 
0  —  {<?i,  •  •  •  j  9n—i}  be  specified  implicitly  in  association  with  the  first 
n  primes  such  that  7  represents  the  average  prime  gap. 

Figure  (5)  is  a  plot  of  Dp(104,  h)  vs.  6  for  all  even  6  from  2  to  104726, 
which  is  the  largest  even  interval  among  the  first  104  primes.  The  points 
are  strongly  concentrated  along  several  distinct  lines  which  converge 
on  the  abscissa  near  5  =  104726,  forming  a  fanned  distribution  in 
the  manner  of  Fig. (3).  We  should  emphasize  that  every  Dx(n,5)  is  a 
singly- valued  function  of  6  for  a  given  n,  and  no  two  points  in  Fig.  (5) 
share  the  same  5.  The  fanned  striations  are  an  aggregate  visual  effect 
that  could  be  produced  only  from  finely-grained  periodicity  in  the  5- 
dependence  of  Dp.  Figure  (6)  is  a  zoomed-in  view  of  Fig. (5)  showing 
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Figure  5.  Distribution  of  intervals  among  the  first  104 
primes. 


non-vanishing  Dp(104,5)  for  6  <  460.  The  quasi-periodic  structure 
observed  in  Fig.  (6)  is  representative  of  the  type  of  structure  observed 
throughout  the  full  range  of  5.  We  have  verified  that  Dp  has  the  same 
basic  form,  with  similarly  well-articulated  periodicity,  for  arbitrary  n 
as  large  as  106,  which  is  the  limit  imposed  by  computational  feasibility. 

Recall  that  the  periodicity  in  Dc  is  a  direct  consequence  of  the  cyclical 
repetition  of  gaps,  which  effects  a  broader  cyclical  repetition  in  the  in¬ 
tervals  spanning  more  than  one  gap.  The  prime  gaps,  however,  are  not 
ordered  in  the  manner  of  repeating  sequence.  It  is  nonetheless  possible 
for  intervals  spanning  many  prime  gaps  to  repeat  quasi- cyclically,  mim¬ 
icking  the  global  properties  of  a  sequence  of  the  form  c.  Specifically, 
suppose  that  there  exists  some  /i  =  n(n)  such  that  the  prime  intervals 
are  globally  subject  to  the  periodic  correlation 


d 


k+/ii 


(33) 
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Figure  6. 

in  approximation  of  the  precise  periodicity  in  Eq.(22).  The  n- dependence 
of  //  has  been  allowed  in  anticipation  of  the  outcome  of  the  following 
analysis.  Let  T  «  represent  the  corresponding  period.  In  order 
for  Eq.(22)  to  be  meaningful  T  =  T{n)  must  correspond  to  a  common 
interval  size  among  the  first  n  primes.  If  the  putative  quasi- cyclical  na¬ 
ture  of  the  prime  intervals  is  strongly  pronounced  then  we  may  expect 
T  ~  Ap.  In  any  case  Eq.(33)  would  imply  that  intervals  of  a  given  size 
5  spanning  l  prime  gaps  occur  with  approximately  the  same  cyclical 
regularity  as  do  intervals  of  size  5  +  T  spanning  l  +  fi  prime  gaps.  We 
would  accordingly  have 

(34)  ggyw±i)  K  gfW) 

n  —  l  —  n  n  —  l 

for  l  <  n  —  2/r,  analogously  to  Eq.(34).  Summation  over  l  produces 


n— 2/z 

IJy.tnJ+Tl  ~  ^ 


n  —  l 
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for  all  8  spanning  on  average  no  more  than  n  —  2 /i  gaps.  Equations 
(33)  through  (35)  could  be  consistent  with  the  PNT  and  the  other 
proven  properties  of  the  primes  while  naturally  effecting  the  periodicity 
observed  qualitatively  in  Dp. 

Ideally  one  should  test  for  the  hypothesized  quasi-cyclical  behaviors 
by  quantifying  the  statistical  strength  of  the  global  correlation  implied 
in  Eq.(33).  The  PNT  guarantees,  however,  that  dj ^  increases  broadly 
with  increasing  j,  for  any  fixed  k—j.  Because  the  putative  oscillations 
associated  with  Eq.  (33)  could  only  be  small  in  comparison  to  the  global 
trend  associated  with  the  PNT,  a  strong  correlation  of  the  form  Eq.(33) 
is  observed  for  nearly  every  meaningful  value  of  p.  A  reliable  direct 
test  of  Eq.(33)  is  therefore  infeasible.  In  contrast,  the  oscillations  in 
Dp  associated  with  Eq.(35)  may  be  large  in  comparison  to  the  broad 
variation  of  Dp  with  respect  to  6.  A  reliable  and  effective  indirect  test 
of  Eq.(33)  may  be  therefore  formulated  by  quantifying  the  strength  of 
the  correlation  implied  in  Eq.(35). 

In  order  to  maintain  consistency  with  Section  2,  the  following  measure 
is  defined  generally  in  terms  of  an  arbitrary  Dx.  We  substitute  for  T 
the  independent  variable  t  and  define  p  ~  i/7  to  be  the  average  number 
of  gaps  spanned  by  an  interval  of  size  t  among  the  first  n  numbers  in 
x.  For  a  given  t  let  8i, ...  ,6m  represent  all  of  the  M  =  M(t)  unique 
values  of  8  in  the  range  (0 ,pn—  2 1\  for  which  Dx(n,  8)  is  non- vanishing. 
With 

(36)  Ui  =  Dx(n,8i  +  t) 


and 


—  n  (  A)  x^Dx\n,Si) 
Vi  =  Dx(n,8i)  -  p2_^ - 

i=i 


(37) 


n  —  l 
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the  standard  coefficient  measuring  the  strength  of  the  correlation  be¬ 
tween  the  left  and  right  sides  of  Eq.  (35)  is 


where  u  and  v  are  the  respective  means  of  iq  and  vt  over  all  <5*.  Finally, 
in  order  to  represent  the  t  for  which  Xx  =  X.x(n,  t.)  is  maximized  for  a 
given  n,  let  rx  =  rx(n)  be  defined  such  that 

(39)  Xx(™,tx)  >  Xx(n,t) 

for  all  t.  If  there  are  more  than  one  t  for  which  yx(n,  t)  has  the  same 
maximal  value  then  let  rx(n)  be  the  smallest  among  them. 

If  the  prime  intervals  repeat  quasi-cyclically  in  the  manner  of  Eq.(33), 
thereby  effecting  periodicity  in  Dp  of  the  form  Eq.(35),  then  yp(n,  t) 
should  exhibit  pronounced  local  maxima  near  unity  for  all  t  equal  to  an 
integer  multiple  of  T,  with  T  typically  producing  the  global  maximum. 
We  therefore  expect  T  ~rp.  If  the  periodicity  is  sufficiently  strong  and 
if  n  is  sufficiently  large  then  we  also  expect  rp«i  Ap. 

Figure  (7)  is  a  plot  of  yp(103,f)  for  even  t  from  6  to  2370.  The  cal¬ 
culations  indicate  a  strong  periodic  correlation  for  all  t  equal  to  an 
integer  multiple  of  6,  and  an  appreciably  less  strong  correlation  for  all 
other  t.  Among  multiples  of  6,  clearly  pronounced  local  maxima  are 
observed  at  all  multiples  of  tp(103)  =  210,  which  produces  the  global 
maximum  xp(103,  210)  ~  0.9914.  We  therefore  conclude  that  the  non¬ 
vanishing  values  of  Zlp(103,5)  and  Dp^(103,<5)  are  mutually  correlated 
in  the  manner  of  Eq.(35)  with  period  T(103)  =  210.  Furthermore  the 
most  common  interval  size  among  the  Erst  103  primes  is  Ap(103)  =210, 
indicating  an  additional  degree  of  consistency  with  the  cyclical  behav¬ 
iors  of  c. 
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Figure  7.  Coefficient  measuring  period-t  correlation  of 
the  form  Eq.(35)  in  Zlp(103,<5). 

We  report  similarly  conclusive  results  for  arbitrary  n  as  large  as  106. 
Figure  (8)  contains  a  plot  of  rp  for  logarithmically  spaced  n  from  102 
to  106.  In  each  case  the  global  maximum  xP(n,  tp)  is  greater  than  0.97, 
and  strongly  pronounced  local  maxima  are  found  at  all  integer  multiples 
of  tp.  Included  in  Fig. (8)  is  a  plot  of  Ap  for  all  n  over  the  same  range. 
Both  Xp  and  Ap  increase  with  n  in  a  remarkably  well  articulated  step¬ 
wise  manner,  remaining  constant  for  a  broad  range  of  n  between  each 
sharp  transition.  Over  the  entire  range  of  n  considered  here  yp  and 
Ap  coincide  all  but  precisely,  differing  only  near  the  transition  points. 
Moreover  all  of  the  observed  values  of  xp  and  Ap  among  the  first  106 
primes  are  of  the  form 

■*(y) 

(40)  y^  =  Y[pa, 

a— 1 

which  is  the  primorial  of  some  real  y  >  2.  The  ‘steps’  in  Fig.  (8) 
correspond  to  =  30,  =  210,  =  2310,  =  30030  and 
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FIGURE  8.  log10  plots  of  the  correlation  period  rp,  in 
heavy  points,  and  the  most  common  intervals  size  Ap,  in 
the  broken  line. 

p7#  =  510510.  Although  yp  is  not  statistically  meaningful  for  n  smaller 
than  approximately  100,  it  is  worth  noting  that  Ap  exhibits  the  same 
step-wise  behavior  for  n  <  100  among  the  primorial  numbers  Pijfz  =  2, 
£>2#  =  6  and  P3#.  We  also  observe  that  each  step-wise  increase  in 
Xp  ~  Ap  from  a  given  p ^  to  pi+ i#  occurs  where  n  is  near  pi+ 1  # . 

The  calculations  presented  here  imply  unequivocally  that  the  non- 
vanishing  values  of  Dp  and  Dp'1  are  correlated  mutually  according  to 
Eq.(35),  with  T  ~rp.  Given  that  Eq.(35)  is  a  specific  and  direct  con¬ 
sequence  of  Eq.(33),  it  is  reliable  to  conclude  that  the  prime  intervals 
repeat  quasi-cyclically  in  the  manner  of  Eq.(33).  The  broad  agreement 
between  rp  and  Ap  reinforces  the  reliability  of  that  conclusion. 

Lastly  it  is  worth  mentioning  that  the  periodic  properties  identified  in 
this  Section  are  associated  with  the  particular  order  of  the  sequence 
of  prime  gaps  and  are  not  determined  simply  by  the  net  distribution 
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G p.  In  order  to  demonstrate  that  point  let  the  sequence  p  =  p(n)  = 
Pi,P2,  ■  ■  ■  ,pn  be  constructed  by  randomly  rearranging  the  sequence  of 
prime  gaps  g-2, .  ■  ■ ,  gn- 1,  with  pi  and  p2  always  equal  to  2  and  3  respec¬ 
tively.  As  a  natural  consequence  of  the  construction  of  p  we  also  have 
Pn  —  Pn-  On  average  a  comparatively  small  number  of  other  pj  happen 
to  be  prime  as  well.  Figure  (9)  is  a  plot  of  Z}p(104,<5)  for  a  particular 
p.  The  present  example  is  representative  of  the  form  of  Dp  for  arbi¬ 
trary  n  and  typical  p.  In  sharp  contrast  to  Dp,  Dp  varies  broadly  with 
6  in  a  simple  linear  manner  and  is  similar  in  form  to  a  distribution 
associated  with  a  sequence  like  r  whose  gaps  are  randomly  selected. 
In  fact  the  linear  trend  of  Z}p(104,<5)  crosses  the  ordinate  near  1900, 
which  is  consistent  with  the  intercept  noj /y  associated  with  Dr  in  a 
case  where  g  is  of  the  form  Eq.(19)  with  uj  —  2  and  7  =  10.47,  where 
10.47  is  approximately  the  average  gap  among  the  first  104  primes.  For 
the  present  purposes,  however,  it  is  sufficient  to  note  that  Dp  generally 
lacks  the  distinctive  periodic  structure  observed  in  Dp. 

4.  Conclusions 

By  examining  the  counting  function  Dp  =  Dp(n,5)  we  have  identi¬ 
fied  several  new  global  properties  of  the  positive  intervals  among  the 
first  n  primes,  for  n  at  least  as  large  as  106.  The  variation  of  Dp 
with  respect  to  5  exhibits  a  strongly  pronounced  periodic  correlation 
of  the  form  Eq.(35),  which  is  the  signature  of  quasi-cyclical  repetition 
in  the  sequences  of  intervals  spanning  fixed  numbers  of  prime  gaps  ala 
Eq.(33).  Over  the  observed  range  of  n  the  characteristic  period  T^rp 
of  the  repetition,  as  inferred  from  the  maximal  correlation  yp(n,  rp), 
almost  universally  coincides  with  the  most  common  interval  Ap.  Such 
a  congruence  is  also  a  signature  of  quasi-cyclical  repetition  of  the  form 
Eq.(33).  Furthermore  for  all  n  considered  here  both  rp  and  Ap  are  al¬ 
ways  primorial  numbers  and  are  typically  the  largest  primorial  number 
no  greater  than  n.  The  regular  manner  in  which  T  increases  with  n 
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Figure  9.  Distribution  of  intervals  in  a  randomized 
analogue  of  the  first  104  primes. 


indicates  a  certain  ‘Russian  dolls’-style  self-similarity  in  the  periodicity 
of  the  prime  intervals,  with  the  characteristics  over  a  certain  range  of 
n  replicated  in  scale  over  consecutively  larger  ranges. 

Although  the  behaviors  reported  here  have  been  validated  only  for  n 
as  large  as  106,  it  is  reasonable  to  expect  that  they  would  be  observed 
for  arbitrary  n.  We  therefore  offer  the  following  conjectures. 

Conjecture  4.1.  For  all  n^>  1  the  positive  intervals  among  the  first 
n  primes  repeat  quasi- cyclically  such  that 

(41)  dj,j+i  ~  dj+^j+i+fj,, 

which  implies 

2/4  (  7“)(0  /  c\  \ 

(42)  Dr(n,S+T)  «  Dv(n,S)  -  ft  j. 
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where  T  —  T(n)  is  the  characteristic  period  and  p  —  /j,(n)  is  the  positive 
integer  nearest  to  the  average  number  of  gaps  spanned  by  an  interval 
of  size  T  among  the  first  n  primes. 

Conjecture  4.2.  For  all  n  >  1  we  have 

(43)  T(n)  «  rp(n)  «  A p(n). 


Conjecture  4.3.  The  most  common  interval  size  A p(n)  among  the  first 
n>  1  primes  is  always  a  primorial  number.  A p(n)  increases  with  n  in 
a  step-wise  manner  from  a  given  ptff  to  Pi+iff,  with  each  transition 
occurring  where  n  is  near  Pi+iff.  We  therefore  have 

(44)  Ap  ~  n. 


Conjecture  4.3  is  stated  without  reference  to  tp  so  that  it  may  be 
considered  independently  of  the  other  Conjectures  and  independently 
of  any  periodic  behaviors  discussed  here. 
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